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Abstract 



Using a method of local transmission matrix, we generalize the well-known 
WKB formula for a barrier penetrability to multi-channel systems. We com- 
pare the WKB penetrability with a solution of the coupled-channels equations, 
and show that the WKB formula works well at energies well below the lowest 
adiabatic barrier. We also discuss the eigen-channel approach to a multi- 
channel tunneling, which may improve the performance of the WKB formula 
near and above the barrier. 
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I. INTRODUCTION 



The coupled-channels approach has been a standard method in describing atomic, molec- 
ular, and nuclear reactions involving internal degrees of freedom [1-5]. In nuclear physics, 
for instance, the coupled-channels method has been successfully applied to heavy-ion fusion 
reactions at energies around the Coulomb barrier in order to discuss the effect of couplings 
between the relative motion of the colliding nuclei and inelastic excitations in the target 
nucleus [1-3]. At these energies, the fusion reaction takes place by a quantum tunneling, 
and the coupled-channels calculations well account for the enhancement of the barrier pen- 
etrability due to the channel couplings. 

A difficulty in the coupled-channels calculations, however, is that it is sometimes not so 
easy to obtain a numerically stable solution with a controlled accuracy. This is particularly 
the case in the presence of closed channels and/or when the coupling strength is strong in 
the classically forbidden region. Several methods have been proposed in order to stabilize 
the numerical solution [4,6-8]. 

In this paper, instead of directly integrating the coupled-channels equations with the 
stabilization techniques, we solve them using the WKB approximation. To this end, we 
employ the method of local transmission matrix, which was originally developed by Brenig 
and Russ in order to stabihze numerical solutions of the coupled-channels equations [4] . We 
solve the equation for the local transmission matrix under the semi-classical assumption, 
and generalize the well-known WKB formula for a barrier penetrability for a single channel 
to coupled-channel systems. Since the penetrability is expressed in a compact form, the 
resultant WKB formula is entirely free from the problem of numerical instability. Moreover, 
the WKB method can be easily applied to systems with a large number of degrees of freedom, 
while obtaining a direct solution of the coupled-channels equations can be computationally 
demanding. Also, the WKB method is useful in gaining a physical intuition for the dynamics 
of multi-channel tunneling. 

The paper is organized as follows. In Sec. II, we set up the coupled-channels equations 
and introduce the local transmission matrix. We derive the semi-classical expression for the 
local transmission matrix, and obtain the WKB formula for a multi-channel penetrability. 
We apply the WKB formula to a three-channel problem, and compare the penetrability 
with the numerical solutions of the coupled-channel equations. In Sec. III. we discuss the 
penetrability at energies near and above the barrier height. The WKB formula which we 
derive works when the multiple reflection of the classical path under the barrier can be ne- 
glected, that is, at energies well below the barrier. We show that the eigen-channel approach 
can provide good prescriptions at higher energies, where the primitive WKB formula breaks 
down. We summarize the paper in Sec. IV. 

II. MULTI-CHANNEL WKB FORMULA 

Our aim in this paper is to derive the WKB formula for penetrability for a one dimen- 
sional potential barrier in the presence of channel couplings. We consider the following 
coupled-channels equations: 
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^Mnno(a;) + [Vnm{.x) + e^5n,m " E] U^no{x) = 0. (l) 

Here, m is the mass of a particle, is the excitation energy for the n-th channel, and E is 
the total energy of the system. Unno i^) is the wave function matrix, where n refers to the 
channel while no specifies the incident channel. Notice that we express the wave functions 
in a matrix form by combining linearly independent solutions of the coupled-channels 
equations, N being the dimension of the coupled equations. For the situation where the 
particle is incident on the barrier from the right hand side, the boundary conditions for 
Unnoix) are given by, 

'^nno('^) ^ T^noC " {x > Oo), (2) 



where A;„ = y2m{E —- e„)/^^ is the wave number for the n-th channel. The inclusive pene- 
trability is then obtained as, 

P = i:^\Tnno\'. (4) 

Let us now introduce the local transmission matrix defined by [4] 

'^(x) = \ [iqixY^u'ix) + u{x)\ , (5) 

where q{x) = [2m{E - W{x))/h'^Y/'^ with Wnm{x) = ^^(x) + e„(5„,„. From Eqs. (2) and 
(3), the asymptotic form of t{x) reads, 

rnm{x) Tnme-'''"'' (x ^ -cx)), (6) 

^ (^n.me"^''"'^ (X^OO). (7) 

Here, we have used the fact qnm(x) — > kn5n,m as x — > ±oo. It is easy to show that the local 
transmission matrix t{x) obeys the equation [4], 

t'{x) = ~q-\x)q'{x)[l + p(x)]t(x) - iqix)T{x), (8) 

where 

p{x) = [iq~'^{x)u'{x) - u{x)] ■ [iq~'^{x)u'{x) + (9) 

is the local reflection matrix [4]. 

The WKB approximation may be obtained by neglecting p{x) in Eq. (8) [9,10], that is, 

t'{x) = -^q-\x)q'{x)T{x) - iq{x)T{x). (10) 

A similar equation has been derived by Van Dijk and Razavy [9,10], but by using the method 
of variable reflection amplitude (see also Ref. [11]). Notice that the asymptotic form of the 
local reflection matrix p{x) is [4] 
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{x — > — oo), 
{x — > oo). 



(11) 

(12) 



Neglecting p{x) in Eq. (8) is thus equivalent to ignoring the reflection, that is reasonable 
in the semi-classical hmit. This, in fact, corresponds to the lowest order of the Bremmer 
expansion [12-16], where the WKB formula is obtained by approximating a smooth potential 
with a series of sharp potential steps ^. 

For a single channel problem, Eq. (10) can be easily integrated to yield 



TiX] 



g(^) g-i qix')dx' 



(13) 



For a coupled-channels problem, however, care must be taken in the integration, since q{x) 
and q'{x) do not commute to each other in general. We attempt to solve Eq. (10) by 
discretizing the coordinate with a mesh spacing of Ax. Replacing the derivative term by a 
simple point difference formula, we obtain 



T{Xn-l) ~ T{Xn) + Ax 



iq{Xn)T{Xn) + -q ^{Xn)q'{Xn)T{Xn] 



1 + \Xn)q'{Xn) 

1 - Axq-'^{xn)q'{xn) 



[l + iq{xn)Ax]T{xn), 



-1/2 



^iq{Xn)^X 



t(x„). 



(14) 
(15) 
(16) 



to the lowest order of Ax. Using q[xn-i) ~ g(a;„)[l — Axq ^{xn)q' {xn)], the first factor in 
Eq. (16) is transformed to be 



[l-Axq-\xr,)q'{xr,)V/^ 



1 



\lq~^Mq{xn-i) ^Jq{xn) 



q{x„). 



(17) 



Approximating 

[q-'^{xn)q{x„-i)]-^^^[q{xn)]~^/^ ~ [q{xn)q~^{xn)q{xn-i)]~^^'^ = l/^q{xn-i), we finally ob- 
tain 



T{Xn-l) 



em-n)^-Jq^x^)r{x^). 



^q{Xn-l) 

Iterating this equation backward from x = oo, we obtain 

1 



(18) 



t(— oo) 



\/9(-oo) 



n 



,iq(xi)Ax 



q{oo). 



(19) 



^The local reflection matrix p{x) satisfies the equation p' = i(qp + pq) — (1 — p)q~^q'{l -\- p)/2 
[4] . One may solve this equation perturbatively assuming that p{x) remains small and finding the 
correction to Eq. (10). 
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Substituting this expression into Eq. (4) together with Eq. (6), the WKB approximation to 
the multi-channel penetrability reads 



n 



n 



jq{xi)Ax 



no 



(20) 



This is the main result in this paper. Notice that the factor kn/kn^ docs not appear in the 
WKB formula, Eq. (20). In practice, one can evaluate Eq. (20) by diagonalizing W{x) at 
each Xi. This yields 



n 



n 



^|m(a;i))e^«-(=^')^^(m(x 



Uq 



(21) 



where \m{x)) is the oigon-vector of the matrix W(x) with the eigen-value of A^(x), and 

qm{x) = \]2m{E - \m{x))/n^. For a one dimensional problem, Eq. (20) is reduced to the 
familiar WKB formula. 



P{E) = exp 



(22) 



where xq and Xi are the inner and the outer turning points, respectively. 

Let us now apply the WKB formula (20) to a three-level problem. We consider the 
following coupling potential: 



W{x) = 




(23) 



with 



Fix) = Poe 



(24) 
(25) 



The parameters are chosen following Ref. [17] to be Vo=100 MeV, -Fo=3 MeV, and s = Sf =3 
fm, which mimic the fusion reaction between two ^^Ni nuclei. The excitation energy e and 
the mass m are taken to be 2 MeV and 29mN, respectively, where mAr=938 MeV is the 
nucleon mass. With these parameters, the three eigen-barriers Xi{x), which are obtained 
by diagonahzing W{x) at each x, have the barrier height of 97.31, 102.0, and 106.7 MeV, 
respectively, while the barrier height for the uncoupled barrier V{x) is 100 MeV (See Fig.l). 

Before we study the three channel problem, let us first examine the validity of the WKB 
approximation for a single channel case to see whether the semi-classical approximation 
works in principle for the parameters which we choose. Figure 2 shows a comparison between 
the WKB penetrabihty for the uncoupled barrier V{x) (the dashed hne) obtained with Eq. 
(20) and the exact solution. It is plotted in the linear and logarithmic scales in the upper 
and the lower panels, respectively. One clearly sees that the WKB approximation indeed 
works well at energies about 2 MeV below the barrier height and lower. 
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As it is well known, the naive WKB approximation breaks down around the barrier. In 
fact, the WKB penetrability is unity at the barrier height, while the exact result is about a 
half. Around these energies, one needs to improve the WKB formula by using the uniform 
approximation in order to take into account the multiple reflection of the classical path under 
the barrier [18,19]. We will discuss this problem more in the next section in connection to 
the penetrability in a coupled-channels system. 

Let us now solve the coupled-channels problem, Eq. (23). We integrate Eq. (21) from 
X = —15 fm to x=15 fm with Ax=0.05 fm. The dashed line in Fig. 3 shows the penetrability 
in the WKB approximation for this problem, which is compared to the exact solution of the 
coupled-channels equations (the solid hne). The figure also shows the penetrability for the 
uncoupled barrier as a comparison (the dotted line). Remarkably, the WKB formula (20) 
reproduces almost perfectly the exact solution at energies well below the lowest adiabatic 
barrier, i.e., 97.31 MeV, as in the single channel problem. 

We notice that the WKB penetrabihty increases much more rapidly than the exact 
penetrability at energies corresponding to the height of each eigenbarrier. This is in close 
analogue to the single channel problem shown in Fig. 2. This behaviour may be expected 
in the eigen-channel approach discussed in Refs. [17,20,21]. We will discuss this point in the 
next section. 



III. PENETRABILITY NEAR AND ABOVE THE BARRIER 

In the previous section, we have shown that the multi-channel WKB formula works 
remarkably well at energies well below the lowest eigenbarrier. Therefore, it can be expected 
that the WKB formula provides a useful framework in discussing, for instance, the role of 
inelastic excitations in the colliding nuclei in nuclear reactions at extremely low energies, 
such as astrophysically relevant reactions, where the standard coupled- channels calculations 
may be difficult to carry out. 

At higher energies, however, we found that the agreement between the primitive WKB 
formula and the exact solution of the coupled-channels equations becomes poor. For a single- 
channel problem, one can cure this problem by using the uniform approximation [18,19]. The 
WKB formula which is valid at all energies is given by, 

PiE) = (26) 



1 + exp [2 dx'^^^{V{x')-E) 



where the turning points Xq and Xi become complex numbers when the energy E is above 
the barrier. It is not straightforward at all, however, to extend the uniform approximation 
to the coupled-channels problem. In this section, we instead present two prescriptions to 
deal with the coupled-channels penetrability at energies near and above the barrier. 



A. Dynamical Norm Method 



The first prescription is closely related to the dynamical norm method developed in Ref. 
[22] . It was argued in Ref. [22] that the penetrability may be expressed as a product of the 
penetrability in the adiabatic limit and a multiplicative factor to it which accounts for the 
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non-adiabatic effect. The latter factor, which was called the dynamical norm factor, was 
evaluated through the imaginary time evolution for an intrinsic degree of freedom with a 
classical path obtained with the adiabatic potential. 

We follow here the same idea as in the dynamical norm method, and re-express Eq. (19) 

as 



(-00) 



q(-oo) 



j~J Qiqo{xi)A.x^i[q{xi}-qo{xi)]Ax 



(27) 



where qo{x) = y2m{E — Xo{x))/?i'^ is the local wave number for the lowest eigen-barrier 
(i.e., the adiabatic barrier), Xo{x). The penetrability P{E) is then given by 



E 



n 



n 



J[q{xi)-qo{xi)]A3 



no 



■ exp 



-2 



xo 




2 (Ao(x') - E) 



(28) 



The second factor on the right hand side (r.h.s.) of this equation is nothing more than the 
WKB penetrability for the adiabatic potential Ao(a:). At this stage, one may replace it by 
the exact penetrability for the adiabatic potential, Pad{E). The first factor on the r.h.s. of 
Eq. (28) expresses the non-adiabatic effect, as in the dynamical norm factor introduced in 
Ref. [22]. Notice that our formula, (28), is in fact an improvement of the dynamical norm 
method in Ref. [22] , since the classical path is not assumed from the beginning in evaluating 
the dynamical norm factor. 

The result of the modified dynamical norm method is shown in Fig. 4 by the dashed line. 
It is evident that the agreement between the WKB approximation and the exact solution 
improves significantly at energies near and slightly above the adiabatic barrier, although the 
method provides essentially the same result as the original WKB formula, (20), at higher 
energies. 



B. Eigen-channel approach 

The second prescription which we discuss is based on the eigen-channel approximation. 
In this approach, the penetrability is expressed as a weighted sum of the penetrability for 
the eigenbarriers [17,20,21], that is, 

P{E)^Y.'>^r^Pn{E), (29) 
n 

where Pn{E) is the penetrability for the eigen-potcntial A„(a;). The weight factors Wn are 
usually estimated by assuming that the matrix W{x) is independent of x through the 
interaction range [17,21], while the coordinate dependence of W{x) is properly taken into 
account in calculating the penetrability, P^i^E). One often takes the barrier position of 
the uncoupled barrier, X^, in order to estimate the weight factors [21]. This leads to (see 
Eq.(21)) 

Wn^Wn{X^)^\{n{X^)\no)\''. (30) 
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This procedure is indeed exact when the intrinsic degree of freedom has a degenerate spec- 
trum [20,23-25], since the matrix W{x) can be diagonalized independent of x. When the 
intrinsic states have a finite excitation energy, however, the unitary matrix which diagonal- 
izes W{x) exphcitly depends on x, and the results may depend strongly on the position 
where the weight factors are evaluated. Also the weight factors possess some energy de- 
pendence in general. In Ref. [20], we have explicitly shown for a two channel problem that 
the optimum weight factors are considerably different from those estimated at the barrier 
position, although their energy dependence appears to be weak. A satisfactory procedure 
to determine the weight factors has not yet been found so far when the excitation energy is 
finite. 

In Fig. 3, we have shown that the WKB penetrability approaches to a constant value at 
the barrier height of each eigen barrier. Assuming that the weight factors are independent 
of energy, one can exploit this fact to determine a consistent value of the weight factors in 
the eigen-channel approach. For instance, at the barrier height of the lowest eigenbarrier, 
E — Bq, assuming that the contribution from the higher barriers is negligible, Eq. (29) 
suggests 

P(So) ~ WoPq{Bo) ~ Wo, (31) 

in the primitive WKB approximation (i.e., without the uniform approximation). Therefore, 
if one evaluates Eq. (20) at E — Eq, it directly provides the weight factor for the lowest 
eigenbarrier. One can repeat this procedure N — 1 times to determine the weight factors 
Wo, Wi, ■ ■ ■ , wn_2- suppose that the weight factors for the k lowest eigen-barriers have been 
determined. The weight factor for the {k + l)-th eigenbarrier is then estimated as 

fc-i 

Wk = PwKBiBk) - (32) 

1=0 

where B^ is the barrier height of the {k + l)-th eigenbarrier, and Pwkb is the penetrability 
in the WKB approximation, (20). Here, we have used the fact Pi{Bf;) = 1 for i < A; in the 
primitive WKB approximation, and assumed Pi{Bk) — for i > k. The weight factor for 
the highest eigenbarrier AAr_i(a;) is evaluated as 

Af-2 

Wn-1 = 1- (33) 
1=0 

in order to ensure the unitarity. 

We apply this prescription to the three channel problem discussed in the previous sec- 
tion. The weight factors are evaluated to be 0.5914, 0.3543, and 0.0543 for the lowest, the 
second lowest, and the highest eigen barriers, respectively. The result of the eigen channel 
approximation with the weight factors thus estimated is denoted by the filled circles in Fig. 
4. The result is indistinguishable from the exact solution of the coupled-channels equations 
for all the energy region shown in the figure. We thus conclude that the multi-channel WKB 
formula which we derive in this paper provides a consistent way to determine the weight 
factors in the eigen-channel approach, and it provides a useful and simple prescription to 
compute the penetrability in the presence of channel couplings at energies from well below 
to well above the potential barrier, as long as the weight factors are slowly varying functions 
of energy. 
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IV. SUMMARY 



We have extended the well known WKB formula for barrier penetrability to systems 
with intrinsic degrees of freedom. Applying the formula to a three channel problem, we 
have explicitly demonstrated that the WKB formula reproduces very nicely the exact solu- 
tion of coupled- channels equations at energies well below the lowest eigenbarrier, i.e., the 
adiabatic barrier. The WKB formula which we derived is applicable even to systems with 
a large number of degrees of freedom, where the standard coupled- channels calculation can- 
not be performed due to a computational limitation. Our method may therefore provide a 
useful framework to discuss, e.g., a quantum scattering problem in the presence of coupling 
to a heat bath [26,27]. The WKB formula is also useful when one discusses the channel 
coupling effect on the tunneling rate at deep subbarrier energies, especially in the presence 
of closed channels, since the direct integration of the coupled-channels equations may suffer 
from a numerical instability. Such interesting problems include heavy-ion fusion reactions 
at extremely low energies [28], electron screening effects in nuclear astrophysical reactions 
[29,30], and nuclear structure effects in astrophysical fusion reactions [31,32]. 

The WKB formula which we derived neglects the effect of multiple reflection of a classical 
path under the potential barrier. Such primitive formula breaks down at energies near and 
above the adiabatic barrier, as is well known. We discussed two prescription to cure this 
problem. One is the dynamical norm method, where the WKB formula is re-expressed as 
a product of the penetrability for the adiabatic barrier and a multiplicative factor which 
accounts for the non- adiabatic effect. By replacing the adiabatic penatrability in the WKB 
approximation by the exact one, we have shown that this prescription improves the result at 
energies near and slightly above the adiabatic barrier. The second prescription is the eigen- 
channel approximation, where the penetrability is given as a weighted sum of the eigen 
penetrability. By applying the WKB formula at energies corresponding to the barrier height 
of each eigen barrier, we have presented a consistent procedure to determine the weight 
factors. Wc have shown that this prescription works well at energies from well below to well 
above the barrier, as long as the energy dependence of the weight factors is negligible. This 
method would provide a useful way to simplify the coupled-channels calculations in realistic 
systems. 

For a single channel problem, the primitive WKB formula can be improved by using the 
uniform approximation [18,19], where the resultant WKB formula is applicable at energies 
both below and above the barrier. It will be an interesting, but very challenging, problem 
to extend it to a multi-channel problem. For this purpose, it will be extremely helpful 
to construct a solvable coupled-channels model. A work towards this direction is now in 
progress [33]. 
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FIG. 1. The eigenpotentials for the three channel problem obtained by diagonalizing the cou- 
pling matrix W{x) at each x (the solid lines). The uncoupled barrier V{x) is shown by the dotted 
line for a comparison. 
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FIG. 2. The barrier penetrability for the uncoupled barrier V{x) as a function of energy E in 
the linear (the upper panel) and the logarithmic (the lower panel) scales. The solid and the dashed 
lines are the exact solution and the WKB approximation, respectively. 
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FIG. 3. Same as fig. 2, but for the three-channel problem. The dotted line shows the penetra- 
bility in the no coupling limit. 
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FIG. 4. Comparison of the barrier penetrability calculated by several methods. The solid line 

shows the exact solution of the coupled-channels equations. The dashed line is obtained with the 
modified adiabatic method which takes the non-adiabatic correction into account, while the filled 
circles are results of the eigen-channel approximation. The upper and the lower panels shows the 
penetrability in the linear and the logarithmic scales, respectively. 
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